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ABSTRACT 

When a star comes within a critical distance to a supermassive black hole (SMBH), immense tidal forces can 
remove a significant fraction of the star's mass, resulting in a stream of debris that falls back onto the black hole 
and powers a luminous flare. In this paper, we perform two hydrodynamical simulations of the disruption of a 
main- sequence star by a SMBH to characterize the evolution of the debris stream after a tidal disruption. We 
demonstrate that this debris stream is confined by self-gravity in the two directions perpendicular to the original 
direction of the star's travel, restricting its width and height to be only a factor of a few larger than its original 
size. As a consequence, the stream has a negligible surface area and makes almost no contribution to either the 
continuum or line emission. We additionally find that the debris stream is strongly compressed and partially 
virialized when it returns to pericenter, resulting in an elliptical structure which grows at a speed somewhat 
less than the star's original orbital speed. We propose that any observed emission lines are not the result of 
photoionization in the unbound debris, but are produced in the region above and below the forming elliptical 
accretion disk, the same region in which the broad-line region is expected to be produced in steadily- accreting 
active galactic nuclei. As each line within a broad-line region is observationally linked to a particular location 
in the accretion disk, we suggest that the absence of a line indicates that the accretion disk does not yet extend 
to the distance required to produce that line. This model can be used to understand the spectral properties of 
the tidal disruption event (TDE) PSl-lOjh, for which Hell lines are observed, but the B aimer series and Hel are 
not. Using a maximum likelihood analysis, we go on to show that a partial disruption of a main-sequence star 
of near- solar composition can reproduce this event, and that the canonical model of an accretion disk truncated 
at twice the tidal radius cannot simultaneously explain the event's luminosity evolution and near-constant color. 
Instead, we propose that some fraction of material that returns to pericenter ensheaths the underlying accretion 
disk, resulting in the majority of the light being converted into a blackbody of near-constant temperature. While 
the formalism developed here is specifically targeted to explain the features of PSl-lOjh, the generalized model 
we present is intended to be the foundational work for classifying future TDEs as they are discovered. 
Subject headings: Keywords here 
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1. INTRODUCTION 

The tidal disruption of a star by a supermassive black hole 
(SMBH) splits the star into either two or three ballistically dis- 
tinct masses. In the event of a full disruption, the star is split 
into two pieces of nearly-equal mass. One half of the star be- 
comes bound to the black hole after the encounter, and contin- 
ues along elliptical trajectories with pericenter distances equal 
to the star's original pericenter distance. The other half of the 
star gains orbital energy in the encounter, and is placed on 
hyperbolic trajectories. For a partial disruption, a third mass 
in the form of a surviving stellar core emerges from the en- 
counter, with the absolute value of its orbital energy compa- 
rable to its own binding energy (Faber et al. 2005; Guillochon 
et al. 2011; MacLeod et al. 2012; Liu et al. 2013). 

Determining the fates of these pieces of the star are crit- 
ical in determining the appearance of the flare that results 
from the immense gravitational energy that will be released 
by the accretion disk that eventually forms. Previously, it has 
been assumed that the unbound material, which was thought 
to be a wide "fan," was the primary contributor to the broad 
emission lines that are produced as the result of a tidal dis- 
ruption (Strubbe & Quataert 2009; Kasen & Ramirez-Ruiz 
2010; Clausen & Eracleous 2011). For the tidal disruption 
event (TDE) PSl-lOjh, it was assumed that hydrogen, which 
is ejected to large distances within the wide debris fan gener- 
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ated by the disruption, can recombine more quickly than the 
rate at which it is ionized by the central source (Gezari et al. 
2012, hereafter G12). This would ensure that the vast major- 
ity of the hydrogen is neutral, and thus any ionizing radiation 
incident upon the fan would produce an emission feature. The 
absence of any hydrogen emission features was used to derive 
an upper limit on the amount of hydrogen present, implying 
that helium is five times more common than hydrogen by mass 
with the disrupted star. 

In this paper, we present three-dimensional hydrodynami- 
cal simulations that show that the assumption that this debris 
fan intercepts a significant fraction of the light is incorrect. 
As noted by Kochanek (1994), the width of the stream of un- 
bound material is still controlled by the stream's self-gravity 
in the transverse direction, restricting its width to only be a 
fraction of the star's original pericenter distance. Through nu- 
merical simulations of fully-disruptive encounters with mass 
ratios q = M^/Mq = 10-^ and 10^, we verify that the transverse 
containment of the stream's width does indeed occur (Figure 
1). As a result, the stream only grows in the radial direc- 
tion, and thus the total volume and surface area increase only 
slightly more steeply than Vp. The surface area of this struc- 
ture is not significant enough to produce hydrogen emission 
lines, even for the disruption of a main-sequence (MS) star 
composed largely of hydrogen. 

But while we find that the area of the unbound debris has 
been vastly overestimated, we also find that the area occupied 
by the accretion disk formed from the bound material has been 
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vastly underestimated. Our numerical simulations confirm the 
prediction that material that returns to pericenter is ballisti- 
cally launched to very large distances from the black hole, 
hundreds of times r^. Additionally, we find that significant 
dissipation occurs when this material returns to pericenter. As 
the debris stream quickly virializes at pericenter and the den- 
sity of the material is significantly reduced as compared to 
the star's original density, self-gravity is suppressed even in 
the transverse direction. As a result, a fan structure is formed 
once material returns pericenter. But as this material belongs 
to the fraction of the original star that is strongly bound to the 
black hole, the radial extent of this material grows at a rate 
that is significantly smaller than the unbound fraction. 

As the region in which Ha is produced in steady AGN is 
on the order of a few light days from the black hole (Peter- 
son 2006), we show that it is very unlikely that the debris 
ejected by the disruption has traveled the distance necessary 
to produce an Ha line for PSl-lOjh. Through comparison 
with the processes responsible for producing the broad line 
regions (BLRs) of steadily- accreting AGN, we predict that the 
helium lines that are observed in PSl-lOjh are produced much 
closer to the black hole (Bentz et al. 2010), and the debris has 
sufficient time to reach this distance by the time the first spec- 
trum was observed. Motivated by the results of our hydro- 
dynamical simulations, we model the accretion disk structure 
and use a Markov-Chain Monte Carlo (MCMC) procedure to 
determine the combinations of parameters with the highest- 
likelihood, and we find that the highest-likelihood models in- 
volve the disruption of a low-mass main-sequence star with 
mass ^O.SMq by a Mh = lO^M© black hole. 

In Section 2 we describe our method for running hydrody- 
namical simulations to characterize the behavior of the debris 
stream after a disruption, and describe the maximum like- 
lihood analysis we employed to estimate the parameters of 
PSl-lOjh. In Section 3 we present a physical interpretation 
of the results of our hydrodynamical simulations. Bearing 
these results in mind, we develop our generalized model for 
the time-dependent, broadband light that would accompany 
the disruption of a star in Section 5. We then apply this model 
specifically to PSl-lOjh in Section 6. Finally, we review ad- 
ditional evidence as to why the disruption of a helium-rich 
star is unlikely to have produced PSl-lOjh in the first place, 
and look towards the future when TDEs will be regularly ob- 
served. 

2. METHOD 
2.1. Hydrodynamical Simulations 

The black hole at the center of our own galaxy is estimated 
to be :^ 4 X 1O^M0 (Ghez et al. 2008), and is one of the small- 
est known massive black holes (Giiltekin et al. 2009). As the 
mass of an average main-sequence star is ~ 0. IM© (Kroupa 
et al. 1993), the majority of stellar tidal disruptions will have 
q> 10^. For such disruptions, the timescale of return of the 
most bound debris is on the order of days to weeks (Rees 
1988), with the peak fallback rate occurring approximately 
one month after the time of the disruption (Evans & Kochanek 
1989; Lodato et al. 2009; Guillochon & Ramirez-Ruiz 2013). 

For hydrodynamical simulations of TDEs, the main limit- 
ing factor is the sound-crossing time of the star, which for 
a solar mass star is approximately one hour. Given an ini- 
tial stellar model that occupies 100^ grid cells, each hydrody- 
namical time- step translates to only one minute of physical 
time. Thus, the simulation of the tidal disruption of a so- 




FiG. 1. — Snapshots from a tidal disruption simulation with M* =Mq, 
Mh = 10^ Mq, and (3 = 1.8, as compared to a simple model of a tidally- 
confined debris stream with self-gravity, where we assume that the width of 
the stream scales as r^/^ (Kochanek 1994). The left panel shows a superposi- 
tion of the debris stream at different times, with the longest stream depicting 
the time when the most bound material returns to pericenter at f = 42 days. 
The color along the stream indicates whether it is bound or unbound from 
the SMBH, with magenta corresponding to bound and cyan corresponding to 
unbound. The green line shows the original path of the star, and the green 
circles show the locations of the surviving core corresponding to the eight 
snapshots shown on the right-hand side of the figure. In each of the right- 
hand panels, the simple model of the debris stream is shown atop the results 
from the simulation. 



lar mass star by a lO^M© black hole that includes the time 
at which the fallback rate is at a maximum requires '^10^ 
time-steps. Additionally, the debris stream resulting from 
the disruption must be fully resolved in both length and in 
width. As the stream is self-gravitating (as described in Sec- 
tion 3.1), it mains a very narrow profile, with the aspect ratio 
of the stream when the first material returns to pericenter be- 
ing ^^^^(^peak/^p)^^^ ~ 10^, where ^peak IS the time where M 
reaches a maximum, and t^ is the pericenter crossing time. 
If the number of grid cells across the stream is forced to be 
at least 20, which is necessary to satisfy the Truelove et al. 
(1997) criteria, then 10^ grid cells would be required to be 
evolved for 10^ time-steps. 

This means that a complete simulation of the full problem 
within a single simulation is very computationally expensive. 
Instead, we run two separate simulations that are each well- 
equipped to describe the behavior of the debris stream at two 
different epochs. To determine the fate of the debris liberated 
from a star during a tidal disruption, we used two similar sim- 
ulation setups, differing only in the mass ratio q. The first 
simulation sets ^ = 10^ and solely focuses on the evolution of 
the debris stream as it expands away from pericenter after the 
star's initial encounter with the black hole. Because of the 
computational expense, the return of the debris to pericenter 
is not followed in this simulation. The second simulation sets 
q= 10^, and follows the return of the debris to pericenter well 
beyond the peak in the accretion rate Mpeak- In these encoun- 
ters, the peak accretion rate is realized only one day after peri- 
center, and we follow the evolution of the returning debris for 
a total of 5 X 10^ seconds (about one week). Our hydrodynam- 
ical simulations were performed in a module written for the 
FLASH adaptive mesh refinement code, the details of which 
can be found in Guillochon et al. (2009, 201 1) and Guillochon 
& Ramirez-Ruiz (2013). 

The initial conditions of the simulation are similar to those 
presented in Guillochon & Ramirez-Ruiz (2013), with the 
poly tropic V that describes the star's structure being set to 5/3, 
and the impact parameter j3 = r^/ri being set to 2. The star 
is placed on a parabolic trajectory at an initial position that 
is several times further than rt, and is initially resolved by 50 
grid cells across its diameter. As realistic equations of state 
are only sensical at the full-scale of the problem, the hydro- 
dynamics of the gas are treated using a simple adiabatic poly- 
trope P (X p^, where 7 is the adiabatic index. The code uti- 
lizes the adaptive-mesh functionality of the FLASH software 
in different ways for the two simulations. In the ^ = 10^ sim- 
ulation, regions which are less than 10~^ times as dense as the 
current peak density are derefined, but maximum refinement 
is maintained within 4rp at all times. For the ^ = 10^ simula- 
tion, each refinement level is assigned to a single decade in p, 
using the star's original central density pc as a baseline, with 
the exception of the first refinement threshold which is set to 
p = 5x 10" Vc 

As both the timescales and the length scales of a ^ = 10^ 
disruption are different from those of the more typical ^ = 10^ 
disruption, care must be taken when interpreting the results 
from these simulations and attempting to scale them up to 
what would be realized for larger mass ratios. As we will 
describe in Section 3.2, the dissipation processes that are ob- 
served in the scaled-down simulation are analogous to other 
dissipation mechanism that operate for larger values of q, 
given the proper scaling. 



2.2. Fitting TDE Observations 

For fitting our models for tidal disruptions to observed 
events, we have developed the code TDEFit, which performs 
a maximum likelihood analysis using an affine-invariant 
MCMC (Goodman & Weare 2010; Foreman-Mackey et al. 
2013). The code is written in Fortran and utilizes the parallel 
variant of the algorithm presented in Foreman-Mackey et al.. 
We have designed the software to be flexible in the model 
parameters it accepts as inputs, any free parameter (either dis- 
crete or continuous) can be included in the parameter space 
exploration by simply listing it and its range of acceptable 
values within a parameter file. In the same way, both trivial 
and non-trivial priors can be specified at runtime for single or 
combinations of input parameters. 

As the solutions can often be multi-modal, with small re- 
gions of acceptable parameter space separated by large voids 
of poor parameter space, it can sometimes be difficult to find 
the deepest global minimum using the vanilla affine-invariant 
algorithm. We have improved upon this algorithm by per- 
forming simulated annealing (SA, Press et al. 1986) on a frac- 
tion F of the walkers every N timesteps during a "bake-in" 
period, where both F and N are adjustable. Each walker that 
is selected to anneal is used to seed an amoeba whose points 
are randomly drawn a small distance away from the original 
walker, these "frackers" then run through a full SA cycle in 
which the temperature is gradually reduced until they are un- 
able to improve upon their local solution. 

This enables the depths of local minima to be found more 
quickly, and in tests we have found that this improves the 
time of convergence to the global solution by orders of magni- 
tude. Additionally, we anneal the ensemble of walkers them- 
selves during the bake-in period, using the temperature sched- 
ule proposed in Hou et al. (2012), and periodically compare 
the scores of walkers to the best so far, removing those that 
fall below a pre-determined threshold that depends on the cur- 
rent annealing temperature. After the bake-in period, the al- 
gorithm reverts to the vanilla affine-invariant MCMC and run 
for several autocorrelation times, ensuring that detailed bal- 
ance is maintained. 

As inputs to this method, we use the full functional forms of 
the fallback rate (parameterized as Msim) presented in Guillo- 
chon & Ramirez-Ruiz (2013) (see their Figure 5) for 7 = 4/3 
and 7 = 5/3 polytropes. We assume that as the mass ratio 
^ ::^ 1, the dependence of M on Mh, M* , and R^ is self-similar. 



M: 



= </'^*,0-*,0 



Rj^^MmW). 
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where 10^Mh,6 = Mh, MqM^^q = M*, /?©/?*,© = /?* and Mint 
is an interpolation of Msim, defined below. We can eliminate 
R^ in this expression by using known mass-radius relation- 
ships (e.g. Tout et al. 1996 for MS stars or Nauenberg 1972 
for white dwarfs). For MS stars, we presume that all objects 
with M* < 0. IMq have the radius of a 0. IM* star. With these 
relations, M is solely a function of Mh, M*, and (3. 

As the simulations of Guillochon & Ramirez-Ruiz are only 
run for specific values of /3, we determine intermediate (3 so- 
lutions by rescaling neighboring simulations in /3- space to the 
same scaled time variable x = (^- Wn)/(Wx- Wn), where ^min 
and ^max are the minimum and maximum times for each Msim 
curve, and then interpolating linearly between the two solu- 
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tions, 



[l3\=rmn{Be(3s,^\B>p} 



(2) 
(3) 



\P]-IP\ 



[M,,^(\p],x)-M,,m([P\.x)], (4) 



where /3sim is the set of all /3 for which a simulation is avail- 
able, and where [/3J and [/3] return the values of /3sim that 
bracket p. We find this preserves the overall shape of the 
M curves well for values of P for which a simulation is not 
available, as long as the sampling in /3sim is sufficiently dense 
to capture the overall trends. 

The objective function used within TDEFit when compar- 
ing our models to the data is the maximum likelihood func- 
tion. 
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where n is the number of datapoints, Vobs,/ and Vmod,/ are re- 
spectively the AB magnitude at the /th datapoint for the ob- 
servation and the model, aobs,/ is the measurement error asso- 
ciated with the /th datapoint, and (jy is the intrinsic variability 
of the source, assumed to be a constant that scales with black 
hole mass (see Section 5.1). 

3. HYDRODYNAMICS OF POST-DISRUPTION DEBRIS 
3.1. Debris stream with self- gravity 

Determining the fate of the various pieces of the star after 
a disruptive encounter is critical in determining the appear- 
ance of the flare that results from the immense gravitational 
energy that will be released by the accretion disk that eventu- 
ally forms. Previously, it has been assumed that self-gravity 
of the disrupted star is unimportant, and therefore the spread 
in energy imparted to the debris at pericenter leads to a spread 
in angle as well as semi-major axis (Strubbe & Quataert 2009; 
Kasen & Ramirez-Ruiz 2010). Under this assumption, the un- 
bound debris is a homologously expanding structure, which 
occupies a constant solid angle and whose volume increases 
proportional to Vp. 

However, Kochanek (1994) showed that the stream can be 
in fact gravitationally confined in the transverse direction by 
self-gravity and forms a very thin structure (Figure 1), with 
a width A and height H that scale as r^/^ for 7 = 5/3, where 
r = r/rt. For general 7, 



H A^-^ (X constant 



(6) 



where A is the mass per unit length (Ostriker 1964), which 
we define to be A = M^/2R^ ait = t^^. Assuming that A ex r~^ 
(Kochanek 1994), 

H = RJ^^, (7) 

where we recover H ex r^/^ for 7 = 5/3. 

In the right eight panels of Figure 1 , we superimpose the 
results of our hydrodynamical simulations of the disruption 
of a star by a black hole with mass ratio ^ = 10^ with this 
simple prescription. We find excellent agreement between the 
prediction of Kochanek and our results over the time period 
in which we ran the simulation. The thinness of the stream 



is also noted in other hydrodynamical simulations in which 
the mass ratio is large (Rosswog et al. 2009; Hayasaki et al. 
2012). 
The surface area of this structure for 7 = 5/3 is 
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where r^ ^ hP^^^it /to) is the distance to which the most un- 
bound r naterial ha s traveled (Strubbe & Quataert 2009), and 
tfi = t:^R%/GM^ is the star's free-fall time. At the peak time 
of ~ 100 days, PSl-lOjh emits '^10^^ erg s~^ of radiation with 
an effective temperature of a few 10^ K, implying a photo- 
sphere size of ~10^^ cm with area '^10^ AU^. By contrast, 
the area occupied by the stream is only comparable to this 
value when t ^ 10^% ^ 10 yr for q = 10^ and P = I. Kasen 
& Ramirez-Ruiz (2010) calculated that this component would 
contribute at most 10"^^ ergs s~^ of luminosity for the disrup- 
tion of a solar mass star. However, we believe that this repre- 
sents an upper limit as the self-gravity of the stream was not 
included in that work, resulting in an artificially fast rate of 
recombination. 

Because the evolution of the stream is adiabatic, but not in- 
compressible, the stream is resistant to gravitational collapse 
in both the radial direction perpendicular to the stream, and in 
the axial direction along the cylinder. Collapse can only occur 
in the radial direction when 7 < 1, as it becomes energetically 
favorable to collapse radially (McKee & Ostriker 2007), and 
can only occur axially for 7 > 2, where the fastest growing 
mode has a non-zero wavelength and leads to fragmentation 
(Figure 4 of Ostriker 1964; Lee & Ramirez-Ruiz 2007). 

In a thin stream, the tidal force applied by the black hole 
results in the density p scaling as r~^ (Kochanek 1994). As 
the distance r oc t^^^ when ^ ^ 00 for parabolic orbits, this 
implies that p oc t~^. For cylinders, the time of free-fall % 
is proportional to p~^^^, the same as it is for spherical col- 
lapse (Chandrasekhar 1961), and thus tff oc t. Therefore, a 
segment of the stream within which % ever becomes greater 
than t will not experience recollapse at any future time, as the 
two timescales differ from one another only by a multiplica- 
tive constant. It is also evident that self-gravitating cylinders 
can be bound whereas a self-bound sphere will not, as the 
Jeans length progressively decreases in size for structures that 
are initially confined in fewer directions (Larson 1985). This 
implies that the cylindrical configuration may not remain self- 
bound if sufficient energy is injected into the star via a par- 
ticularly deep encounter in which the core itself is violently 
shocked, which only occurs for /3 > 3 (Kobayashi et al. 2004; 
Guillochon et al. 2009; Rosswog et al. 2009), approximately 
10% of disruption events. Additional energy can also be in- 
jected by nuclear burning (Carter & Luminet 1982; Rosswog 
et al. 2008), again only for events in which P is significantly 
larger than 1. For most values of P, the amount of energy in- 
jected into the star at pericenter is insufficient counteract grav- 
itational confinement in the transverse direction at r = 0, and 
thus from the timescale argument given above the unbound 
stream would forever be confined. 

Given this result, and given the computation burden of re- 
solving a structure with such a large aspect ratio for the num- 
ber of dynamical timescales necessary for the material to be- 



gin accreting onto the black hole, we presume that the gravi- 
tational confinement continues to hold at larger distances than 
we are capable of resolving, and show how the profile of the 
debris stream would appear if the simulation were followed to 
the point of the material's return to pericenter in the left panel 
of Figure 1. 

By contrast, the bound material travels a much shorter dis- 
tance from the black hole before turning around. For the ma- 
terial that remains bound to the hole, it has been assumed that 
the material circularizes quickly after returning to pericenter, 
resulting in an accretion disk with an outer radius equal to 2rp, 
where r^ is the pericenter distance (Cannizzo et al. 1990; Ul- 
mer 1999; Gezari et al. 2009; Lodato & Rossi 2010; Strubbe 
& Quataert 2011). This is actually a vast underestimate of 
the distance to which the debris travels, which can be found 
via Kepler's third law for the orbital period of a body and di- 
viding by two to get the half-period, and then solving for the 
semi-major axis a, 






1/3 



(9) 



where we have made the assumption that r^ c^ 2a, appropri- 
ate for the highly elliptical orbits of the bound material (The 
most-bound material has eccentricity e = l-2q~^^^ = 0.98 for 
^ = 10^). As this material is initially confined by its own grav- 
ity, the return of the stream to pericenter mimics a huge /3 
encounter, which as we explain in the following section can 
yield a impressive compression ratio. 

3.2. Dissipative effects within the nozzle 

There are a combination of potentially active mechanisms 
that can provide the required dissipation for any given event, 
with each mechanism dominating for particular combinations 
of Mh and /3. To quantify the effect of each of these mech- 
anisms, we define the ratio V = (dE/dt)(T /E), where T is 
the orbital period. V represents the fraction of gravitational 
binding energy that is extracted per orbit, with V = 1 indicat- 
ing a mechanism that fully converts kinetic to internal energy 
within a single orbit. To have M and L trace one another over 
the duration of a flare, V must have a value > ^m/^peak, where 
^peak is the time at which the accretion rate peaks and tj^ is 
the time at which the most bound debris returns to pericenter. 
For partially disruptive encounters, the ratio between these 
two times is ~ 3, but then increases as /3^ for deep encoun- 
ters in which tj^ varies more quickly than ^peak (Guillochon & 
Ramirez-Ruiz 2013). 

In the following sections we provide a brief description of 
the dissipation mechanisms that are expected to operate in a 
TDE. Only the first mechanism (hydrodynamical dissipation) 
is present within our calculations, as we do not include mag- 
netic fields or the effects of a curved space time. Regardless 
of the origin of the dissipation, we expect that the dissipation 
observed in our simulations is likely to be quite analogous to 
the other dissipative process that operate. 

3.2.1. Hydrodynamical dissipation 

As described in Section 3.1, self-gravity within the debris 
stream sets its width and height to be equal to R^r^^^, and thus 
when the stream crosses the original tidal radius, its height is 
approximately equal to the size of the original star. If the re- 
turn of the material to pericenter behaved in the same way 
as the original encounter, the maximum collapse velocity v^ 



would be equal to the sound speed at rt multiplied by (3, yield- 
ing a dissipation per orbit V = q~^^^/P, equal to 2% for ^ = 10^ 
and /3 = 2 (Carter & Luminet 1983; Stone et al. 2012). In our 
hydrodynamical simulations f or ^ = 10^, we find that ~ 10% 
of the total kinetic energy of the debris is dissipated upon its 
return to pericenter through strong compression at the nozzle 
(Figure 2). This is a factor of a few larger than the expected 
dissipation. 

However, as the star has been stretched tremendously, the 
sound speed within the stream has dropped by a significant 
factor, meaning that the distance from the black hole at which 
the stream's sound-crossing time is comparable to the orbital 
time (i.e., where the tidal and pressure forces are approxi- 
mately in balance) is not the star's original tidal radius, but 
is instead somewhat further away. 

For two points that are separated by a distance dr within 
the original star, their new distance d/ upon returning to peri- 
center is related to the difference in binding energy between 
them, which remains constant after the encounter and is equal 
to 

dr~rr,~ Rl 



-q 



1/3 



(10) 



As angular momentum is approximately conserved, the two 
points will cross pericenter at the same location they origi- 
nally crossed pericenter, but at two different times t separated 
by a time dt owing to their different orbital energies. Assum- 
ing that the star originally had approached on a parabolic or- 
bit, and that all the bound debris are on highly elliptical orbits, 
the distance from the black hole is 



-GM^f 



1/3 



and thus 



dr' _ fAGM^V'^ 
~dt~ \ 3t J 



(11) 



(12) 



If we set ^ = to be the time when the first point re-crosses 
pericenter, the difference in time is simply the difference in 
orbital period, and thus we can use Kepler's third law to derive 
dE/dt. Using Equations 10 and 12, we can use the chain rule 
to determine dr /dr. 



dr' 
dr 



dE dt dr' 
dr dE dt 



2/3 



4/3 



<y2/3 



(13) 



For q= 10^, this implies that the material has been stretched 
by a factor > 10^ at all times beyond the time the most-bound 
material begins accreting, and for ^ = 10^ this factor is > 10"^. 
The change in volume is then given by the change in cross- 
section of the stream multiplied by the change in length given 
by Equation 13, 

dV dr' , 



dV dr 



-y-y-l 



(14) 



where we have used Equation 7 to estimate the width and 
height of the stream. 

The density of the stream p as it returns to pericenter can 
be approximated by assuming that dM/dr = A, although in 
reality this distribution can be determined more exactly from 
the numerical determination of M by a change of variables 
from E to r. Under this assumption, the change in density 
is simply related to the change in volume alone. As the tidal 
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Fig. 2. — Columndensity contours of the debris resulting from the disruption of a M* =Mq star by a Mh = 10^ Mq black hole. The column density shown in all 
panels is scaled to the value that would be expected for a disruption by a Mh = 10^ Mq black hole, which is 10~^-^ smaller than what is expected for Mh = 10^ M© . 
The six mini-panels in the upper right show the evolution of the column density in the x_y-plane with time, with the upper left mini-panel showing the column 
density at the time of return of the most bound material tm, and lower left mini-panel showing the column density at the time of peak accretion fpeak- The three 
large panels show the column density as viewed from the x-y, x-z, and y-z planes at t = fpeak- 

where 7 is the adiabatic index of the fluid. Atr = r^^^, the ratio 
of the sound speed within the stream Cs,s to the star's original 
sound speed Cs,* is 



radius is proportional to p^/^, the ratio of the effective tidal 
radius of the stream rt,s to rt is then 






dV 



Stt 
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2 7-1 

3 47-5 






(15) 



where we have substituted rt^s/h for r, and where we have pre- 
sumed that the time until the stream reaches pericenter from 
rt,s is small compared to the time since disruption. 

Under the assumption that the stream expands adiabatically 
and its pressure is governed by ideal gas pressure, this results 
in a reduction in the sound speed c^ = ^/dP/dp ex y(i~^V2^ 
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fdV'\ 
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(16) 



Analogous to the original star, the maximum collapse ve- 
locity of the stream vx,s is equal to the sound speed at rt,s 
multiplied by the stream's impact parameter (3^ = /3rt,s/rt, 
v± = /5s<^s,s- As the majority of the dissipation comes through 
the conversion of the kinetic energy of the vertical collapse via 



shocks, the fractional change in the specific internal energy is 



Vh^ 



2 2 
s ^s,s 
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2 (7-1X37-5) 
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4 (7-1X37-5) 
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M 



-2/3 



(17) 



As rt (X p^/^, and c^ oc p ^/^ for 7 = 5/3, Equations 15 and 
16 are inverses of one another in the adiabatic case, 



'^ = '^ = eoM'J''(^ 



(18) 



and thus 17 simplifies to Vh = Pq~^^^ , identical to the amount 
of dissipation experienced by the original star. One key dif- 
ference exists between the original encounter and the stream's 
return to pericenter: While the original encounter may only 
result in the partial shock-heating of the star, even for rela- 
tively deep P (Guillochon et al. 2009), the fact that the col- 
lapse of the stream is highly supersonic (/3s ~ 60/3) guaran- 
tees that shock-heating will occur upon the material's return 
to pericenter. 

For our ^ = 10^ simulation, the amount of dissipation ex- 
pected per orbit predicted by Equation 17 is 4 x 10~^, and for 
our ^ = 10^ simulation the expected dissipation would only be 
4 X 10""^. Thus, the conversion of kinetic energy to internal en- 
ergy via shocks at the nozzle point is inefficient for all but the 
lowest mass ratios and/or the largest impact parameters, and 
would be incapable of circularizing material on a timescale 
that is shorter than the peak timescale of PSl-lOjh (Figure 3). 
This suggests that a viscous mechanism that involves an un- 
resolved hydrodynamical instability, or a mechanism that is 
beyond pure hydrodynamics, is responsible for the circular- 
ization of the material for this event. 

An additional complication that is not addressed here is re- 
combination. As the stream expands, its internal temperature 
drops below the point at which hydrogen and helium begin to 
recombine, flooring its temperature to ^10"^ K (Roos 1992; 
Kochanek 1994). This implies that the ratio between the ini- 
tial and final sound speeds is somewhat smaller than when 
assuming adiabaticity holds to arbitrarily-low stream densi- 
ties, and depends on the initial temperature of the fluid, which 
is ~ 10^ K in the outer layers of the Sun, but ~10^ K in its 
core. This also causes the stream to expand somewhat due 
to the release of latent heat. However, as the material returns 
to pericenter, the compression of the material will reionize it. 
Given these complications, it is unclear if this process would 
lead to more or less dissipation at the nozzle. 

3.2.2. Dissipation through General Relativistic Precession 

For orbits in which the pericenter is comparable to the 
Schwarzschild radius rg, the orbital trajectory begins to de- 
viate from elliptical due to precession induced by the curved 
space-time. The precession time in the inner part of the disk 
is (Valsecchi et al. 2012): 



7gr = 



t) c2 1-^2' 



(19) 



where T and e are respectively the period and eccentricity of 
the stream. As the debris resulting from a tidal disruption has 
a range of pericenter distances (rp ±7?*), there is a gradient 




Fig. 3. — Fraction of binding energy dissipated at f = fpeak for three mech- 
anisms that may contribute to the circularization of material after a tidal 
disruption, with light green corresponding to hydrodynamical shocks at the 
nozzle point, light blue corresponding to dissipation through GR precession 
(presuming 7 = 5/3), and white corresponding to the MRI mechanism. For 
each mechanism, three contours of V are shown, with solid corresponding 
to 100%, dashed corresponding to 10%, and dotted corresponding to 1%. 
If all three mechanisms operate, the shaded blue regions represent zones in 
which V adopts the values specified by the unions of the regions enclosed 
by the three sets of contours, with the lightest/darkest corresponding to the 
least/most dissipation. In order for L and M to trace one another closely, 
V ~ 100%. For reference, our two hydrodynamical simulations are shown 
by the cyan triangle and the red hexagon, and the highest-likelihood fit re- 
turned by our maximum likelihood analysis (Section 6.3) is shown by the 
magenta square. 

in precession times of the returning debris. This precession 
causes the orbits to cross one another, dissipating energy (Er- 
acleous et al. 1995). When compared to the standard a viscos- 
ity prescription, the timescale of this precession is comparable 
to the viscous time. 



W = 10"'''^677'/'(l+^) 



yr, 



(20) 



where 75 = lO^T is the local disk temperature. In a tidal dis- 
ruption, the most-bound material is also the material with the 
shortest precession time, and it is this timescale that sets the 
overall rate of dissipation. By setting T and e in Equation 19 
to t^ = ^/q/2l3~^tf{ and ^^ = 1-2/3^"^/^, the period and ec- 
centricity of the most-bound material, and by assuming that 
precession through an angle 27r would lead to complete dis- 
sipation, the dissipation due to relativistic precession for the 
material that corresponds to the peak in the accretion rate is 



V( 



GR,peak - 



_ ^peak 



2^x2/33^ 



M' 



2/3 



1-^2 



(21) 



In general, the period of the most-bound material tends to 
smaller values for larger q and /3, resulting in more dissipa- 
tion, except in the case that rp and R^ are comparable (Figure 
3, cyan curves). If none of the other dissipation mechanisms 
are effective, this means that disruptions by massive black 
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holes, in which r^ and r^ are closer to one another in value, 

would be the only cases in which M and L follow one another 
closely. 

As we will show in Section 6, the highest-likelihood mod- 
els of PSl-lOjh seem to be consistent with a relativistic en- 
counter with the central black hole. If no other dissipative 
mechanisms are able to operate effectively in tidal disruption 
accretion disks, the possibility exists that only those events in 
which Tp and r^ are comparable will have L and M trace one 
another closely. 

3.2.3. Super-Keplerian, Compressive MRI 

The magnetic field in the bound stellar debris is likely to be 
amplified by compression and by strong shearing within the 
nozzle region. The effects of magnetic shearing, the magneto- 
rotational instability (MRI; Balbus & Hawley 1998), is ex- 
pected to lead to the rapid exponential growth of the magnetic 
field with a characteristic timescale of order the rotational pe- 
riod. This instability has been routinely studied in the context 
of accretion disks and we argue here that it is likely to oper- 
ate within the nozzle although its exact character is difficult 
to constrain given that the boundary conditions are constantly 
changing and the material is never in steady-state. 

The MRI is present in a weakly magnetized, rotating fluid 
wherever 

-^<0. (22) 

dmr 

The ensuing growth of the field is exponential with a charac- 
teristic time scale given by ^mri = A'K\dVt / d\nr\~^ (Balbus & 
Hawley 1998). For a (super-)Keplerian angular velocity dis- 
tribution n (X r~^/^, this gives ^mri = (4/3)Q~^. Exponential 
growth of the field on the timescale Q~^ by the MRI is likely to 
dominate over other amplification process such as field com- 
pression within the same characteristic time. While a vari- 
ety of accretion efficiencies are reported in numerical realiza- 
tions of magnetically-driven accretion disks, which depend on 
the geometry, dimensionality of the simulation, and included 
physics, essentially all models find that the strength of the 
magnetic field amplifies to the point that it is capable of con- 
verting fluid motion into internal energy. From global simu- 
lations of the MRI, the build-up of the magnetic field strength 
is confirmed to be exponential, resulting in a dissipation rate 
that is a constant multiple of the orbital period. In Stone et al. 
(1996), this constant is found to be 3. Once the magnetic field 
strength is saturated, the resulting angular momentum trans- 
port will be governed by the turbulence and is therefore ex- 
pected to take place over longer timescales. A simple estimate 
of the saturation field can be obtained by equating the charac- 
teristic mode scale, ~ VAidft/dlnr)'^, where va is the Alfven 
velocity, to the shearing length scale, ~ dr/dlnft, such that 
^sat ^ (47rp)^^^ftr. This saturation field is achieved after tur- 
bulence is fully developed, which in numerical simulations 
takes about a few tens of rotations following the initial expo- 
nential growth (Hawley et al. 1996; Stone et al. 1996). 

For the Sun, the initial interior magnetic field energy at the 
base of the convective zone Eb^o ^ 10~^^£^g (Miesch & Toomre 
2009), although larger initial fields are possible in general 
(Durney et al. 1993). As the tidal forces stretch the star into 
a long stream, the volume of the fluid increases by a factor 
(3^ (Equation 15) prior to returning to pericenter, reducing the 
magnetic field strength further. However, when the stream re- 
turns to pericenter, it experiences a dramatic decrease in vol- 



ume by a factor /3s (Luminet & Carter 1986). Assuming 

the frozen flux approximation, the new magnetic energy den- 
sity is 

Eb = Eb^oP?^ (23) 
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(24) 



For 7 = 5/3, the dependence on (3^ disappears, i.e. the mag- 
netic field strength upon return to pericenter is identical to the 
star's initial interior field. Assuming that the magnetic field 
within the debris has strength Eb relative to the local gravita- 
tional binding energy Eg upon returning to the nozzle, Vmri 
adopts a simple form (Figure 3, white curves). 



J^MRI,peak - TT ^^P 



^peak 

3L" 



(25) 



3.3. Is the Debris Disk Dissipative Enough? 



In order for the emergent luminosity L to follow the feed- 
ing rate M closely, the dissipation must be effective enough 
such that material returning to pericenter can circularize on 
a timescale tc that is at most the time since disruption ^d- In 
Figure 3 we show the three sources of dissipation that we es- 
timated above. We find that while significant dissipation is 
expected for large (3 encounters, or for encounters in which 
Tp ~ Tg (as is the case for our best-fitting model for PSl-lOjh), 
that there are many combinations of (3 and q that may not 
have the required dissipation necessary to ensure the direct 
mapping in time of M to L. In our ^ = 10^ hydrodynamical 
simulation, we found somewhat more dissipation than what is 
expected from a simple analytical calculation, but the resolu- 
tion at which we resolved the compression at pericenter was 
only marginally sufficient to resolve the strong shocks that 
form there. 

One potential resolution to this issue is the adiabatic index 
of the fluid 7, which in the above calculations we have as- 
sumed = 5/3, although the real equation of state within the 
stream is likely softer due to the influence of recombination. 
With a softer equation of state, the cancelations that occur for 
7 = 5/3 and eliminate the dependence on (3^ for the hydro- 
dynamical (Equation 17) and MRI (Equation 25) dissipation 
mechanisms would no longer apply, yielding both increased 
compression and magnetic field strengths, and thus additional 
dissipation. 

While the initial dissipation of the stream may indeed come 
as the combination of the three previously described mech- 
anisms, it is likely that the mechanism responsible for the 
accretion onto the black hole once the material has been as- 
sembled into a disk is the MRI mechanism, as is suspected 
for steadily-accreting AGN. Given the computational chal- 
lenge of simultaneously resolving the nozzle region and the 
full debris stream, it is clear that local high-resolution mag- 
netohydrodynamic simulations are required to determine the 
true dissipation rate V at the nozzle. 

4. THE RELATIONSHIP BETWEEN STEADILY-ACCRETING AGN 
AND TDE DEBRIS DISKS 

Within the debris structure formed from a tidal disruption, 
the same mechanisms that operate in steady-state AGN may 
continue to operate. There are a number of differences be- 
tween the structure of a debris disk resulting from tidal dis- 
ruption and the structure of steadily-accreting AGN, but we 



will argue that similar processes are responsible for the ap- 
pearance of both structures. In this section, we will make 
continued reference to the highest- likelihood model of PSl- 
lOjh, which is determined in Section 5. 

4.1. The Conversion of Mass to Light 

For steadily-accreting AGN, energy is thought to be re- 
leased by the viscous MRI process at all radii. The amount 
of energy available at a particular radius depends on the lo- 
cal gravitational potential, and thus the vast majority of the 
energy emitted by accreting black holes is produced within a 
few times the Schwarzschild radius r^. The temperature pro- 
file that results from this release in energy within the accre- 
tion disk is given by the well-known expression presented in 
Shakura & Sunyaev (1973), and scales as r"^/"^, resulting in a 
sum of blackbodies with a continuum slope FyOzu^^^ (Pringle 
& Rees 1972; Gaskell 2008). AGN are divided into two pri- 
mary categories: Compton-thick AGN, which are obscured by 
^10^"^ cm^ of column and thus making them optically thick to 
Compton scattering (Treister et al. 2009), and thermal AGN, 
which have column densities significantly less than this value, 
enabling the black hole's emission to be directly observed. 
However, all thermal AGN show an excess in the blue known 
as the "big blue bump" (Shields 1978; Lawrence 2012), and 
the slopes of their continuum iv ex y~^ (Gaskell 2009). This 
is more consistent with the notion that the light emitted from 
the central parts of the disk is intercepted by intervening gas 
before it is observed, with nearly one hundred percent of the 
light emitted by the disk being reprocessed in this way. This 
implies that a significant fraction of the mass that may eventu- 
ally be accreted by the black hole is suspended some distance 
above the disk plane, where it can intercept a large fraction of 
the outgoing light. 

For the accretion structure that forms from the debris of a 
tidal disruption, the dissipation at the nozzle point provides a 
means for lifting material above and below the orbital plane, 
resulting in a sheath of material that surrounds the debris and 
is very optically thick for certain lines of sight (Figure 2). 
However, as the spread in energy at the nozzle point does 
not completely virialize the flow, the resulting distribution of 
matter is flattened, allowing the central regions of the accre- 
tion disk to be visible through material that is close to the 
Compton-thick limit. The time- series presented in the upper 
six panels of Figure 2 show that despite the continually ac- 
tive dissipation process at the nozzle-point, the region directly 
above the central parts of the accretion disk remain relatively 
evacuated of gas before, during, and after the time of peak 
accretion. For the toy simulation presented here, the optical 
depth to Thomson scattering directly above the black hole and 
perpendicular to the orbital plane of the debris is ~ 1, depend- 
ing on the electron fraction of debris (assumed to be pure hy- 
drogen in Figure 2). For more massive black holes, the debris 
is spread over a larger volume, as the tidal radius grows as 
M^^ , and thus N ex M^^ . If the dissipation rate were the same 
independent of black hole mass, it would be expected that the 
disruption of stars by more massive black holes would yield 
more lines of sight for which r ~ 1 . 

4.2. Source of Broad Emission Lines 

Broad line emission is visible in many AGN, being thought 
to be produced by gas above and below the disk plane at dis- 
tances on the order of light days away from the black hole. 
For other AGN, this region is not directly observable, which 



has been attributed to a torus at large radii that can obscure the 
broad line region for lines of sight that run within a few tens of 
degrees of the disk plane (i.e. the AGN unification model, An- 
tonucci 1993). The emission lines produced within this region 
have been successfully used to measure black hole masses 
(Peterson & Wandel 1999; Peterson et al. 2004) based on mea- 
surements of the time lag in the response of line luminosity to 
variations in the output of the central engine. It is still debated 
whether this material is in the form of an optically-thick disk 
wind (Trump et al. 2011) or optically-thick clouds (Celotti & 
Rees 1999), but in either case the material that constitutes the 
BLR is likely bound to the black hole (Proga et al. 2008; Pan- 
coast et al. 2012). 

In a steady-accreting AGN, material accretes from very 
large distances (> 10^ rg), and the emission from this region 
is often manifest as an IR bump in Type II AGN (Koratkar 
& Blaes 1999). At such distances, the ionizing flux originat- 
ing from the black hole is not sufficient to maintain a large 
ion fraction within the disk's emitting layer. The closer one 
gets to the central black hole, the greater the incident flux of 
ionizing radiation on the BLR wind/clouds that generate the 
observed emission lines. The fraction of atoms in an excited 
state X+ relative to the state directly below it Xq is approxi- 
mately (Frank et al. 2002) 



nx^ ^ ajyion) 6(Xo) 

^Xo aB(Xo)/l^on47rr2^e 



(26) 



where 2(Xo) is the flux in photons capable of ionizing the 
lower state. This expression shows that that as the distance 
from the central engine increases, the number of atoms in the 
high state decreases, assuming that the electron density n^ de- 
creases with radius more slowly than r~^ (as Q(xr~^), and also 
shows that species with larger ionization potentials will have 
less atoms in the high state than species with smaller ioniza- 
tion potentials. This leads to a hierarchy of ions in the disk, 
with those with the highest ionization potential being predom- 
inant in the inner regions of the disk. In a steadily- accreting 
AGN, the flux in ionizing photons is large enough to fully ion- 
ize iron (as evidenced by the existence of Fe K lines, Fabian 
et al. 2000), and given that atoms at large radii are mostly 
neutral, all ionic species of all elements exist at some distance 
from the central engine. Reverberation mapping supports this 
basic photoionization picture, as Rblr ^ L^^^ (Bentz et al. 
2010, 2013). In particular, the optical wave band hosts sev- 
eral lines from the B aimer series of hydrogen and lines from 
both singly-ionized and neutral helium (Bentz et al. 2010). 

This wide range of scales is in stark contrast to the de- 
bris disk formed as the result of a tidal disruption, which we 
schematically illustrate in Figure 4. Rather than material spi- 
raling in from parsec scales, material is instead ejected from 
the nozzle point, which lies at the star's original point of clos- 
est approach, and typically has scales on the order of a few 
AU. As a result, the debris disk forms from the inside out. 
The ratio of line strengths in a TDE is dependent upon the 
number of atoms in the photosphere that are in the particu- 
lar ionization state associated with each line. For PSl-lOjh, 
the lack of an Ha emission line was interpreted by G12 as 
being attributed to a lack of hydrogen atoms. However, the 
B aimer series requires neutral hydrogen to be present in suf- 
ficient quantities to produce a line in excess of the continuum 
emission. As shown in our ^ = 10^ simulation, material is 
ejected from the nozzle point at approximately the escape ve- 
locity, with the fastest moving material traversing a distance 
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Fig. 4. — Schematic figure from our q = 10^ simulation demonstrating the geometry of the debris resulting from a tidal disruption Sitt-td= 4.3 x 10^ s, shown 
from the top (top image) and the side (bottom image). The three-dimensional isodensity contours are colored according their temperature, with red being hot 
and blue being cold. Super-imposed on these contours are a line of arrows showing the path of the circularizing debris as return to the black hole (black disk), 
inside which a surviving stellar core may reside (white disk). The dot-dashed and double-dot-dashed lines respectively show the regions interior to which helium 
is doubly-ionized and hydrogen/helium are singly ionized within the BLR. The region in which BLR clouds may form is super-imposed using the gray contours, 
although we note that the BLR may instead be in the form of a diffuse wind. 



^t[(^~^d)/^p]^ • This sets an upper limit on the radial extent of 
the disk. Therefore, the lack of an observed emission feature 
may simply be the result of the disk not being large enough to 
host the region required for that particular feature's produc- 
tion. 

The specifics as to which particular radii contribute the 
most to the emission strength of each line is complicated to 
determine, and requires a more-through treatment of the ion- 
ization state of the gas as a function of radius, which depends 
on the geometry of the structure, and the distribution of den- 
sity and temperature as functions of height and radius. In 
a tidal disruption, the matter distribution that ensheathes the 
black hole is established quickly, forming a steady- state struc- 
ture that is supported by a combination of gas pressure and 
angular momentum (Loeb & Ulmer 1997). Accretion then 
proceeds through the midplane, in which the majority of light 
is generated within a few r^ at X-ray temperatures. These 
photons are intercepted by the ensheathing material at higher 
latitudes. Korista & Goad (2004) determined the equivalent 
widths of various lines as functions of volume density and 
ionizing flux, which is not expected to vary much as a func- 
tion of column density for 10^^ < A^ < lO^^ cm^ (Ruff 2012). 
In Figure 5 we show a series of profiles corresponding to our 
highest-likelihood model for PSl-lOjh over the range of times 
at which spectra were taken of this event. To calculate the 
density distribution ^(H) as a function of r, we once again use 
the chain rule. 



X(H) dMdE f2a\ 
47r m^r^ dE da \r^ J 



3/2 



(27) 



where X(H) is the mass fraction of hydrogen, and presuming 
that the radial distribution of mass is determined by the dis- 
tribution of mass with semi-major axis a, dM/da, set at the 



time of disruption. Likewise, dM/dr is directly proportional 
to dMjda, with a scaling factor equal to the ratio of time spent 
at apocenter versus pericenter, dM/dr c^ dM/da{2a/r^)^/^, 
where we have assumed that 1 - ^ ^ and thus the apocenter 
distance r^c:i2a. 

As a strong dissipation mechanism likely operates at the 
nozzle point, and this dissipation mechanism is likely to be 
as dissipative as the commonly invoked MRI mechanism, it 
stands to reason that the vertical structure of the debris disk 
formed through the circularization process is similar to that 
of a steadily-accreting AGN. Therefore, we would expect that 
the BLR associated with such structures should be similar to 
the BLR produced by steadily-accreting black holes. Under 
this assumption, we can use Equation 27 to approximate the 
number density of hydrogen as a function of radius and to 
determine the equivalent width of various emission lines us- 
ing the models that have been generated for steadily accreting 
AGN (Korista & Goad 2004). Figure 5 shows the density pro- 
files calculated from Equation 27 as a function of time for our 
highest-likelihood model of PSl-lOjh, with the purple curve 
corresponding to the time of the first acquired spectrum of 
PSl-lOjh at -22 days, and the red curve corresponding to the 
last recorded spectrum at +358 days. 

From Figure 5, it is clear that the equivalent width of Hell 
A4686 is significantly larger than that of Ha, H/3, and Hel 
A5876, all three of which are not observed in PSl-lOjh. The 
figure does suggest that hydrogen and/or singly-ionized he- 
lium emission lines may appear at later times when the ioniz- 
ing flux has decreased, although this may not ever be observ- 
able in PSl-lOjh where the flux originating from the TDE has 
already dropped below that of the host galaxy. 

In making this plot, we have made some assumptions that 
actually would lead to a decrease in the strength of the un- 
observed lines if we performed a more-detailed calculation. 
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Fig. 5 . — Contours of the log of the equivalent width of four emission lines 
as a function of hydrogen-ionizing flux <I>(H) and hydrogen number density 
w(H), where the black dashed and solid contours correspond to 0.1 and 1 
decade, respectively (Adapted from Figure 1 of Korista & Goad 2004), with 
the smallest contour corresponding to 1 A of equivalent width. The colored 
triangle within each panel indicates the peak equivalent width for each line. 
The rainbow-colored curves show the profiles of the debris resulting from 
the tidal disruption that corresponds to the highest-likelihood fit of PSl-lOjh 
(Section 6.3), with the range of curves being shown corresponding to the 
range between the first and last spectrum taken for the event, with purple 
being -22 and red being -1-358 days from peak. The black dotted curve shows 
the conditions at Tq over the full event duration. Note that Ha, H/3, and Hel 
A5876 would potentially be observable if additional spectra were collected at 
later times. 

Firstly, the models of Korista & Goad (2004) presume that a 
full annulus of material exists at each radius; this is not the 
case in an elliptical accretion disks where the inner annuli are 
closer to full circles than outer annuli (Eracleous et al. 1995). 
In fact, it is unlikely that the outer material can circularize at 
all, given that there is significantly less angular momentum 
in the disk than the angular momentum required to support 
a circular orbit at the distance at which these lines would be 
produced (at r = 10^^ cm, ~30 times more angular momen- 
tum would be required to form a circular orbit than what is 
available at rp). Secondly, we have made the assumption that 
the material that does the reprocessing remains at the distance 
determined by the energy distribution set at the time of dis- 
ruption at all times (a la Loeb & Ulmer 1997), when in real- 
ity the entire debris structure will shrink onto the black hole 
due to dissipation at pericenter. It is possible that this shrink- 
age of the debris could prevent emission features arising from 
species with lower ionization potentials from ever being ob- 
served. 

Radiation pressure (which we ignore in this work) may act 
to push some fraction of the material outwards, which in prin- 
ciple could produce low-energy emission features (Strubbe 
& Quataert 2009). However, our highest-likelihood models 
predict a peak accretion rate that is sub-Eddington, and thus 
only a small fraction of the accreted matter is expected to be 
driven to large distances via radiation. It is unclear whether 
the amount of mass in this component would be dense enough 
to produce these features, as the recombination time may be 



too long. 

5. A GENERALIZED MODEL FOR THE OBSERVATIONAL 
SIGNATURES OF TDES 

As emphasized in the previous sections, there are many un- 
certainties relating to how the material circularized when it 
returns to pericenter, how this returning material radiates its 
energy when it falls deeper into the black hole's potential well, 
and in the ionization state of the gas within the debris super- 
structure. Using the code TDEFit, developed for this paper, 
we construct a generalized model of the resultant emission 
from TDEs. In this section, we describe the results of running 
this fitting procedure, and how PSl-lOjh specifically allows 
us to evaluate some of the other models that have been pro- 
posed for modeling TDEs. 

5.1. Model description and free parameters 

Our generalized model for matching TDEs is one in which 
an accretion disk forms by the disruption of a star of mass M* 
by a black hole of mass Mh with impact parameter /3 and off- 
set time ^off = ^-^d, where ^o is the time the first datapoint was 
collected. This disk spreads both inwards and outwards from 
Tp, and is ensheathed by a diffuse layer of material that inter- 
cepts some fraction of the light. The disk itself is bounded by 
an inner radius n and outer radius To, with ri assumed to be set 
by the viscous evolution of the material, and Tq being set by 
the ballistic ejection of material as it leaves the nozzle region, 
which scales as To = r^itjt^^l^ , where ^m is the time of return 
of the most-bound material. The fraction of the full annu- 
lus di that is covered by the disk varies as a function of time, 
with 6>f = when r = To, and = 27r when t = tyi^c (r), assuming 
its spread in the azimuthal direction is controlled by the local 
value of the viscosity. The model is shown pictographically in 
the left panel of Figure 6, with the aforementioned size scales 
as functions of time being shown in the right panel of the same 
figure. 

The source of this viscosity may be similar to the source 
of viscosity at the nozzle point (see Section 3.2), or it could 
be the result of the stream-stream collision that occurs when 
material reaches apocenter (Kochanek 1994; Kim et al. 1999; 
Ramirez-Ruiz & Rosswog 2009). For simplicity, we assume 
that the same viscous process, parameterized by the free pa- 
rameter V, applies in both regions. 

The emergent emission from the disk is calculated using 
the prescription of Done et al. (2012), which largely follows 
the original prescription of Shakura & Sunyaev (1973), but 
amends the no-torque boundary condition to include the ef- 
fects of the black hole's spin, parameterized by the dimen- 
sionless spin parameter ^spm. The model accounts for a shift in 
the emergent disk spectrum arising from variations in opacity, 
resulting in an effective temperature that can be ~2.7 times 
larger than expected from the fiducial SS model. We do not 
permit the luminosity L to exceed the Eddington luminosity 
^Edd = 47rGM^mpc/K,i, and set L = L^dd at times where M ex- 
ceeds this limit. We also include the inclination of the struc- 
ture relative to the observer (/) as a free parameter, where ^ = 
is defined to be edge-on, assuming that both the disk's height 
and the ensheathing layer scale with V in the same way, with 
the emergent emission from both components being reduced 
by a factor V + (1 - V) cos (cj)). 

Note that we assume the color correction is intrinsic to the 
disk emission, and is not the same as the reprocessing that 
occurs due to the diffuse gas that ensheaths the disk and is 
ejected from the nozzle region directly. The photosphere of 
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Fig. 6. — Evolution of the size scales relevant to the appearance of TDEs. The left image shows a three-dimensional schematic of the elliptical accretion disk 
(rainbow-colored surface), the low-density material that ensheaths the disk (cyan surface), and the location of the region interior to which helium is doubly- 
ionized (magenta surface), at three times labeled A, B, and C. The right plot shows results from the highest-likelihood fit of PSl-lOjh, where the solid cyan and 
dashed magenta curves correspond to the two surfaces in the left image, the dash-dotted blue curve corresponds to the distance to which the unbound material 
has traveled, and the dotted orange curve corresponds to the distance to which light travels since the time of the accretion of the most-bound material, denoted by 
the vertical dashed black line. The times to which the images on the left correspond are shown by the labeled vertical black lines. 



this reprocessing layer, whose size is set by a combination of 
the mass distribution and the absorption process responsible 
for intercepting the light, is less constrained. For steadily- 
accreting AGN with thermal emission, the reprocessing layer 
has temperatures of several 10^ K (Koratkar & Blaes 1999; 
Lawrence 2012), and intercepts nearly 100% of the emission 
from the disk, resulting in an effective photosphere size that 
can be hundreds of AU in size for SMBHs accreting at the Ed- 
dington limit. However, as there are many non-thermal AGN 
whose spectra are more representative of bare slim-disk mod- 
els (Walton et al. 2013), it remains unclear how the size of this 
reprocessing zone and its fractional coverage are set. 

In the models of Strubbe & Quataert (2009, 2011), this 
layer is presumed to arise as the result of ejection via a super- 
Eddington wind, and scales with this value when M exceeds 
^Edd- From our hydrodynamical simulations, we find that the 
material that forms the reprocessing layer may be deposited 
by a process that does not require the accretion rate to exceed 
^Edd. but instead depends on the details of how energy is in- 
jected into the material within the nozzle region. The distribu- 
tion of mass in radius resulting from the ejection from the noz- 
zle maps is directly related to dM/dE, although it is modified 
somewhat by the additional spread in energy introduced at the 
nozzle point. However, this spread in energy is local to mass 
that return at a particular time t, and thus the distribution of 
mass with radius after leaving the nozzle point will resemble 
dM/dE with an additional "smear" equal to the spread in en- 
ergy applied at the nozzle. Therefore, we expect that the mass 
distribution with radius follows the general shape of dM/dE, 
and that there will be a density maximum corresponding to 
the orbital period of the material that constitutes the peak of 
the accretion. This peak in density that corresponds to the 
apocenter of the material that determines Mpeak is clearly seen 
in our hydrodynamical simulations (Figure 2). 

Thus, the size of the reprocessing layer is likely to be depen- 
dent on both the instantaneous value of M, which determines 
the amount of ionization radiation produced by the disk and 
the rate of instantaneous mass loss from the nozzle region, 
and on the integrated amount of mass that has been ejected 
from the nozzle region since t = t^. The optical depth r is 



Jo 



dM 
dr 



where k. is an opacity that we leave as a free parameter, but is 
at minimum the Thomson opacity n^ = 0.2(1 -\-X(H)) cm^ g~^ . 
As we described in the derivation of Equation 27, the 
amount of mass at a particular distance r is related to the 
amount of mass at a particular binding energy E, and thus 
we can rewrite Equation 28 in terms of £^, 



Eo 



dM 

Ie' 



-dE. 



(29) 



where Ej^ and Eq are the binding energy of the most bound 
material and the material at apocenter at time t respectively. 

We presume that the fraction of light C intercepted by the 
reprocessing layer scales simply with r. 



C=l- 



(30) 



but enforce the condition that C < c(t-td)/Rph (where Rp^ is 
the size of the photosphere) at all times, otherwise the photon 
diffusion time would be greater than the time since disruption, 
and thus L and M would not be expected to closely trace one 
another. In reality, C should have a wavelength dependence, 
but for the purposes of this work we treat the opacity as being 
"gray," absorbing all frequencies of light equally. 

As the ionization state of the gas (and therefore the opacity) 
depends on the current luminosity, the size of the photosphere 
is expected to vary with time. In general, as the Thomson 
cross-section is significantly smaller than that of bound-free 
transitions, the photosphere scale is likely to correspond to the 
first species is not completely ionized, and in the case of PSl- 
lOjh where Hell emission is observed, we speculate that this 
species is helium (Figure 6, magenta region). If we assume 
that 7?ph (X M^ Tph (X M'^, and L (X M (X ^p^^ph. then the power 
law indices of Rp^ and Tph are simply related. 



2/ + 4m=l. 



(31) 



Jr, 



(28) 



If the opening angle of the reprocessing layer is independent 
of r, the flux in ionizing photons intercepted is constant, im- 
plying 1 = 1/2 and thus m = 0, i.e. Tph is independent of time. 
However, as we find that the geometry may in reality be some- 
what more complicated (Figure 2), we do not assume the in- 
tercepting area necessarily scales as M, and instead leave / 
as a free parameter. For any / ^ 1/2, the temperature of the 



13 



photosphere will evolve with time. We leave / as a free pa- 
rameter and relate m and / through Equation 3 1 . The size of 
the photosphere is then defined to be 



^ph - ^ph,0<^p Tv — 
\MEddJ 



SGM^f 



^peak ^rr 

TT 



1/3 



(32) 
(33) 



where a^ is the semi-major axis of the material that accretes 
ait = ^peak, and 7?ph,o is a dimensionless free parameter. 

The amount of reddening in the host galaxy is also an un- 
known quantity that must be fitted to simultaneously with 
the parameters of the disruption. For extinction in the IR 
through the UV, we adopt the reddening law fits of Cardelli 
et al. (1989), in which the amount of reddening is defined by 
A(A) = Av[a(A) + /7(A)/7?v], in which a(X) and b(X) are fitted 
parameters, Ry is a fitted parameter that ranges between 2 and 
10 (Goobar et al. 2002), and where we take A^h = 1 -8 x lO^Uy 
g cm~^. For the X-rays, we adopt the cross-sections presented 
in Morrison & McCammon (1983). Extinction in the X-rays 
is particularly sensitive to metallicity and temperature (Gnat 
& Ferland 2012), and the uncertainty in the amount expected 
for a particular event is large given the environment of a galac- 
tic center is likely to have super-solar metallicities (Cunha 
et al. 2007) and a wide range of temperatures and densities 
(Quataert 2002; Cuadra et al. 2006; De Colle et al. 2012). 

An advantage of the MCMC method employed here is that 
it permits the inclusion of discrete parameters that can only 
assume particular values. This enables us to simultaneously 
fit multiple physical models, as long as the continuous pa- 
rameters are shared between the models. We include two dis- 
crete free parameters in this work: A^ , which parameterizes 
the type of object that was disrupted, and A^, which param- 
eterizes the polytropic model that is assumed. We include 
two distinct object types, the white dwarf sequence and the 
main sequence. Within each of these sequences, different 
mass ranges are characterized by different polytropic 7; we 
use the M functions derived from our hydrodynamical simula- 
tions (Guillochon & Ramirez-Ruiz 2013) appropriate to each 
mass range. 



, V , V , . / MS : 



M = M5/3(0< 



MS : 
WD : 



O.3<M*/M0 <22 
M,/Mq > 1.0 

M^/Mq < 1.0 
M^/Mq > 22 
M^/Mq < 1.0 



(34) 
(35) 



where some overlap is permitted in the mass range 0.3 < 
M^/Mq < 1.0 to account for the gradual transition between 
fully radiative and fully convective stars in this range. It was 
found that the white dwarf sequence, which is only permits 
very low mass black holes, is excluded to very high con- 
fidence for all of the combinations of parameters that were 
considered, especially when accounting for the measurement 
of the black hole's mass presented in G12, which restricts 
M > 2 X 1O^M0. For simplicity, we exclude discussion of 
the white dwarf channel for the rest of this work. 

In addition to modifying the emergent disk spectrum, ^spin 
also affects the minimum approach distance of a star on a 
parabolic trajectory (i.e. the innermost bound circular orbit) 
^iBCO. and the spread in energy across the star at pericenter 
(Kesden 2012). For simplicity in this work we only consider 



prograde encounters (^spm > 0), and apply first-order correc- 
tion factor to the binding energy. 



E'=ll-Y-f^) 



-1/2 



(36) 



In general, retrograde and/or orbits in which the orbit's incli- 
nation is not equal to the black hole's spin inclination would 
be expected. 

AGN show variability from the radio to the X-ray, with vari- 
ability on the order of a few tenths of a dex being common 
(Webb & Malkan 2000). While the photometric errors are 
small for this event, it is clear that the light curve exhibits 
some intrinsic variability, as may be expected for an accret- 
ing black hole. To model this, we add an additional intrinsic 
spread a^ in quadrature with the observational errors asso- 
ciated with each data point. In addition, the variability has 
been shown to be dependent on the black hole mass (Uttley 
& McHardy 2005). To include this within our modeling, we 
define a function based on the break in the power spectral den- 
sity (PSD) of black hole variability using the results of Kelly 
etal. (2011) 



= 0.0253: 



+0.071 »^-0.19±0.78 



-0.038 



m; 



(37) 



where ^h is the timescale of the break in the PSD, and ( is 
the square root of the variability amplitude measured at the 
break. The quoted error bars are calculated using the fits of 
Kelly etal. 
In total, our fitting procedure includes 15 parameters, 13 of 



which are continuous (M*, Mh, /3, ^off, ^^spm, V 



h ^ph,0, 



7?v, A/h, and (jy), and 2 of which are discrete (A^ and A^). 

5.2. Using the existence/absence of emission lines and their 
properties to constrain TDEs 

In addition to using the quality of the fit of the model light 
curves to the data, we also impose additional constraints de- 
pending on which lines do or do not exist in spectra taken at 
various times (Figure 7). To do this, we measure the emergent 
flux at A = 5100, which is used in steadily- accreting AGN to 
measure the continuum, and compare To to the distance im- 
plied by the relationship between ALa(5 100) and 7?blr, as first 
determined by Wandel et al. (1999) for H/3. Since then, the 
relationship between L and /?blr has been more-accurately 
determined for H/3 (Peterson et al. 2004; Bentz et al. 2013), 
and for several other emission lines (Bentz et al. 2010). For 
all lines, it is found that L ex ^blr' indicating that the source 
of ionizing photons is point-like and that the disk maintains a 
relatively constant scale-height for a wide range of r, which 
gives the natural result that the number of ionizing photons 

Our modification to the likelihood function is simple: If a 
line exists in a spectrum and To < /?blr, or if a line doesn't 
exist and To > ^blr, we reduce the log-likelihood measured 
from the light curve alone >Clc by a factor 



ln/:BLR = ln 



2 \ 2crBLR 



(38) 



where ctblr is the error in the measured L-7?blr relation, 
which we take from Bentz et al. (2010). 

If a line exists, and a velocity for that line has been mea- 
sured, we can use that additional information to constrain the 
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Fig. 7. — Location of BLR in AGN as a function of L at A = 5100 A as compared to the outer radius of the truncated debris disk resulting from a tidal disruption. 
The purple error bars and purple curve show the data and fit of Peterson (2006) for H/3, while the red and blue curves show the offsets to the best-fitting H/3 
relation as measured by Bentz et al. (2010) for Ho; and Hell A4686. The outer radius of the truncated debris disk Vo is shown with the black curve, with the time 
of collection for the four spectra (-22, -1-227, -1-254, and -1-358 restframe days) being denoted by the black points. In all but the last spectrum, in which no lines 
are apparent. Hell A4686 was observed and Ho; was not observed, which implies that Vq must extend beyond the region within which the majority of the flux for 
Hell A4686 is produced, but not extend beyond the region within which most of the Ho; flux is produced. Thus, the four black points must appear between the 
red and the blue curves to satisfy this constraint (shaded in gray). Additionally, the velocity dispersion of the Hell A4686 line (cthcIi) was measured for the first 
spectrum taken at -22 days. The Hell A4686 line is assumed to be located along the standard BLR relation, and the velocity of this line predicted by the model is 
given by matching the luminosity to its associated size on the Hell A4686 curve (white point). As it is unclear whether the motions being observed are Keplerian 
or inflow/outflow, the distance implied by the model velocity is constrained between a circular and a parabolic orbit, <J^h/f^Heii ^ ^Heii < ^GMh/a^^jj (dashed 
vertical blue curves, shaded in blue). 



event further by relating 7?blr to the underlying velocity ex- 
pected at the position. In steadily-accreting AGN, there is 
an ambiguity in the velocity measured for the BLR and its 
location, as it may partially be supported by radiation pres- 
sure (Marconi et al. 2008), we ignore this effect here as it 
only introduces a '^10% correction (Barth et al. 2011). A 
greater uncertainty exists in TDE debris disks in that the un- 
derlying velocity vblr can range from Keplerian (vblr = vk) 
to parabolic (vblr = a/2vk), and therefore we cannot constrain 
the distance implied by vblr better than a factor of >/2. Bear- 
ing this in mind, our reduction to the log-likelihood assumes 
the following functional form. 



ln£v= < 



V<VBLR |(™)' 

Vblr < v < V^vblr 



(39) 



V > Vblr 



f V- V2vblr ] 



We show the results of imposing these constraints for PSl- 
lOjh in Figure 7, where we plot the distance to which material 
has traveled Vq versus the luminosity at 5100 A. The specific 



constraints we have applied are that Hell A4686 must be pro- 
duced, and Hof must not be produced, in the four spectra in 
which the observed light is not dominated by the host galaxy 
(at -22, +227, +254, and +358 restframe days). We find that r^ 
is sufficiently large to produce Hell A4686 in all four of these 
spectra, and that Ha would potentially be observable at later 
times if the host galaxy did not dominate the observed light 
(the light is already subdominant to the host galaxy at +254 
days). H/3, which is produced at smaller distances than Ha 
in steadily-accreting AGN, may potentially be observable in 
late-time spectra, but its wavelength (4861 A) notably over- 
laps with the broad Hell A4686 feature, and is usually a fac- 
tor of several weaker than Ha in most AGN (Vanden Berk 
et al. 2001). We also might expect that H7 and/or Hel A5876 
may appear at later times, as the distances at which these 
lines are produced are only slightly larger than the distance at 
which Hell A4686 is produced (Bentz et al. 2010). As we had 
mentioned in Section 4.2, we may also be overestimating the 
size of reprocessing region if it accretes onto the black hole 
quickly, which would tend to predict the existence of more 
emission features. Given these uncertainties, we do not im- 
pose a constraint on the non-existence of Hel A5876, H/3, or 
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Fig. 8. — Fits of highest Hkelihood of PSl-lOjh for two different models: A bare SS disk (left panel, described in Section 6.2), and a disk with a variably-sized 
reprocessing layer (right panel, described in Section 6.3), with the AB magnitude of the data and the models shown in the top panels, and the difference between 
the highest-likelihood model and the data shown in the bottom panels. The highest-likelihood model found is shown by the solid curves, whereas the 2-a range 
in magnitudes encompassed by the full ensemble of walkers is shown by the shaded bands. The five colors correspond to four filters of Pan-STARRSl system 
(denoted by Pg, Pi, Pr, and Pz), and the NUV band of the GALEX instrument (denoted as GN). 



H7 in this work; we note that their inclusion would likely re- 
strict the size of the accretion disk further, which would tend 
to favor lower-mass black holes. 

We note that the constraints we are imposing do not con- 
sider the specific luminosity of the lines versus the continuum 
(see Section 4.2), which can strongly affect whether a line is 
identified within a collected spectrum. This means we also 
cannot consider in detail the effects of the elliptical accretion 
disk structure resulting from a tidal disruption on the strengths 
of the observed lines, which would preferentially reduce the 
strength of lines originating at large distances as the BLR does 
not occupy a full 27r in azimuth in the outskirts of the debris 
structure (Figure 2). 

6. MODEL FITTING OF PSl-lOjh, A PROTOTYPICAL TIDAL 
DISRUPTION 

6.1. Available Data 

For the fitting procedure, we use all of the available data 
to constrain the event, including four Pan-STARRS bands (Pg, 
Pr, Pi, Pz), the X-ray upper limits from the Chandra space- 
based X-ray telescope (cycle 12), and the spectra taken by 
the Hectospec instrument on the MMT telescope, all of which 
are taken from G12. As the data presented in G12 is already 
corrected for extinction assuming A/h = 7.2x 10^^ cm~^, we 
remove this correction before using the data as an input, as we 
self-consistently determine the extinction in the model fitting 
process. We assume a redshift z = 0.1696 as is determined in 
G12 from template fitting to the host galaxy. 

6.2. Bare Disk Model 

In steadily-accreting AGN disks, the majority of radiation 
produced by the disk is thought to be intercepted by inter- 
vening gas that reprocesses the original emission from the 
disk. In a TDE, this layer may take some time to form, or 
may not form at all, depending on the dissipative processes 
at work. In this case, the light produced as the result of a 
TDE would resemble a bare Shakura-Sunyaev (SS, Shakura & 



Sunyaev 1973) or slim-disk (Abramowicz et al. 1988) model, 
with peak emission that extends well beyond the tens of eV 
that is characteristic of AGN spectra. For this model, we do 
not include the additional constraints imposed by the exis- 
tence/absence of emission lines. 

Fits assuming this "bare" disk model, which are equivalent 
to our generalized model for TDE debris disks sans the repro- 
cessing layer (cyan region of Figure 6), are shown in Figure 
8, where it is immediately obvious that a bare disk model is a 
poor match to PSl-lOjh. In order for the bare disk model to 
closely match the data, the fitting routine settles upon one of 
two non-ideal solutions: An SED with peak that centers about 
the range of wavelengths covered by the observed bands, or 
an SED in which the bands are all within the Ray leigh-j cans 
tail. In the first case, the luminosity L can closely follow M, 
but the color evolves tremendously as the peak of the summed 
blackbody curves shift into/out of the observed bands. In the 
latter case, the ratio of fluxes between the observed bands re- 
mains constant, but L scales as a much weaker power of M, 
L (X M^/4 (Figui-e 9). 

It has previously been assumed that the size of the disk is 
controlled by the angular momentum content of the returning 
material, which is limited to y/2GM^i (see Section 3.1 for 
references). In this case, the disk resembles a bare disk that is 
"truncated" at r = 2rt. We find that these models are an even 
poorer match to the event. 

6.3. Fits to Generalized Model With Reprocessing Layer 

From the previous section, we know that bare disk models 
can either reproduce a constant color, or reproduce a lumi- 
nosity that follows M, but cannot reproduce both behaviors 
simultaneously. This suggests that a secondary process is in- 
volved that reprocesses a large fraction of the light prior to 
reaching the observer. In section 6.3 we suggested that this 
mechanism is the absorption of the soft X-ray photons pro- 
duced primarily at r ~ rg by material deposited at r ~ ^p. 

As can be seen in the right panels of Figures 8 and 9, these 
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Fig. 9. — Spectral energy distribution of the emergent radiation from the highest-hkehhood fits of PSl-lOjh as a function of time for the bare disk (left panel) 
and generalized (right panel) models. The colored curves are spaced equally in ln(f -^d), with red corresponding to early times and blue correspond to late times. 
The vertical black lines show the wavelength of highest transmission for the four Pan-STARRS and GALEX NUV bands. The bare disk models feature a single 
component that corresponds to the accretion disk that peaks at A ~ 10^ A, whereas the generalized model features two components which corresponds to the 
accretion disk, which peaks at A ~ 10^ A, and the single-temperature reprocessing component, which peaks at A ~ 10^ A. 

the host galaxy, but rather the large intrinsic scatter in the Mh- 
L relation, as black holes of mass 10^ < M\^/Mq < 10^ have 
been found for other galaxies of similar magnitude (Graham 
& Scott 2013). The disruption is predicted to have occurred 
42 days prior to the first observation, about 20 days prior to 
what was originally suggested in G12. 

We find that ^spm is only loosely constrained, with the main 
effects of a larger spin being that deeper- /3 encounters would 
be permitted (which are disfavored anyway), and an increase 
in the efficiency of converting mass to light. The inclination (j) 
is highly degenerate with this parameter, and shows a strong 
anti-correlation (i.e. more-slowly spinning black holes tend to 
be more face-on). We find that the preferred models increase 
V to as large a value as possible, and likely this result would 
be altered given a physical model for V that accounts for all 
the various dissipation processes (see Section 3.2). 

The mass-averaged opacity n = 0.5 cm^ g~^ is slightly larger 
than Thomson, consistent with the expectation that the opac- 
ity is dominated by electron scattering with some contribution 
from bound-free absorption by partially-ionized species heav- 
ier than hydrogen (with potentially a large contribution com- 
ing from singly-ionized helium). We find that / = 0.96, which 
is larger than 1 /2 and indicates that the photosphere temper- 
ature does evolve to some degree during the encounter, with 
T (X M~^-^^ (Equation 31). The fact that the highest-likelihood 
models settle on this temperature evolution is a consequence 
of the increase in flux in the Pz band relative to other bands at 
late times. A truly isothermal photosphere could be achieved 
if its location is mediated by the recombination of a single 
species (acting as a thermostat), in which case m = 0. The 
fact that the highest-likelihood models deviate from this value 
indicates that either the distribution of matter in radius and 
height may be non-trivial, or that multiple species may be 
responsible for setting the observed temperature (Lawrence 
2012). The photosphere scale parameter /?ph,o = 0.17 corre- 
sponds to a physical size of 5.4 x 10^"^ cm at peak (approxi- 
mately 50 times larger than pericenter distance rp = 1.1 x 10^^ 
cm), the time evolution of which is shown in the right panel of 
Figure 6. This distance is approximately the distance at which 
helium becomes doubly ionized (Figure 5). 

For the extinction in the host galaxy, we find that a column 
of A^h = 10^^ cm~^ is preferred, with a reddening parameter 
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models provide excellent fits to the data; the parameters asso- 
ciated with the fits of highest likelihood are shown in Table 1 . 
We immediately caution the reader that the reported medians 
of the probability distributions, and the small spread in distri- 
butions of some parameters, should not be taken at face value. 
In our generalized model, which is only a simplified realiza- 
tion of the true structure of the debris, we have made many 
assumptions, and the uncertainty in some of these assump- 
tions is likely to be greater than spread of solutions about the 
highest-likelihood models presented here. That being said, it 
is encouraging that such a simple model with relatively few 
free parameters can provide a reasonable fit to the data, and is 
highly suggestive of the true values of the underlying param- 
eters. 

In the generalized model, we find that the disruption is best 
matched by the nearly-complete disruption (/3 = 0.87) of a low 
mass star (M* = O.56M0) by a Mh = 1.1 x \QPMq black hole. 
This combination of parameters is close to the most common 
sub-Eddington disruption expected (De Colle et al. 2012), but 
predicts the black hole mass is a factor of a few larger than 
the black hole mass suggested by G12. Most of this discrep- 
ancy is likely to arise not from improper template fitting of 
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Fig. 10. — Posterior distributions of Mh, M*, flspin. and ^ph,o for PSl-lOjh. 
Within each panel are the probability P scaled to the maximum probability 
^max, with the red curves showing the posteriors when the BLR constraints 
are not included, and the blue curves showing the posteriors when they are 
included. With the constraints, slightly smaller Mh and slightly larger M* are 
preferred, leading to a slightly larger photosphere scaling parameter /?ph,o to 
compensate. 

Ry = 6.5. This value is somewhat higher than what is typi- 
cally observed within the Milky Way (Ry = 3.1), and is more 
representative of "gray" dust in which all wavelengths are ab- 
sorbed equally. Such values of Ry have been observed outside 
of the Milky Way (see e.g. Falco et al. 1999), and are typical 
of dense molecular clouds (Draine 2003). Another possibility 
that our generalized model simply does not produce enough 
UV photons, necessitating a gray opacity law to compensate. 

Lastly, we find that the model requires ay = 0.034 magni- 
tudes of intrinsic variability, well within the expected range 
given a black hole mass of lO^M© (Equation 37), and surpris- 
ingly small given the potentially chaotic nature of the accre- 
tion process. This strongly suggests that the accretion process 
is smooth and regular, with no major changes in global struc- 
ture over short timescales. 

We find that our highest-likelihood models with and with- 
out the BLR constraints are very similar to one another. In 
Figure 10, we present the posteriors of four parameters that 
showed the most change (Mh, ^*, ^spin, and Rphfi) when the 
BLR constraints are removed, and find that difference in the 
posteriors is on the same order as the scatter about the me- 
dian. This suggests that the timescale, luminosity, and color 
of PSl-lOjh are sufficient to constrain most of the physical 
parameters of an event, whereas the BLR constraints can be 
used as a sanity check to ensure there is no discrepancy be- 
tween the existent BLR emission regions and the observed 
spectra. 

7. DISCUSSION 

7.1. Arguments against the helium star interpretation 

The discovery of a flare with no noticeable hydrogen fea- 
tures certainly hints at the possibility that the disrupted star 
may have been relatively devoid of hydrogen. Aside from the 
hypothesis presented in the previous sections, there are other 



reasons to believe why the helium-rich progenitor scenario 
might be unlikely. 

Firstly, helium-rich stars are rare in the universe. The 
known candidates are SdB/SdO stars (~10^ in the MW, Han 
et al. 2003), helium WDs (--10^ in the MW, Nelemans et al. 
2001), and WR stars (-10^ in the MW, van der Hucht 2001). 
While there is some evidence that the mass function around 
SMBHs is not well-represented by a canonical IMF (Bartko 
et al. 2010), it seems unlikely that the numbers of these stars 
could be increased by the factor of ~ 10^-10^ required to plau- 
sibly explain why the first well-resolved TDE happened to be 
a helium-rich star. 

A second possibility is that the helium-rich star comes as 
the result of the previous interaction of a giant star with the 
SMBH, or potentially through a collision between the giant 
and a more compact stellar object (Davies et al. 1991). How- 
ever, in both of these cases, hydrogen is not completely re- 
moved from the star. In fact, even for deep tidal encounters, 
the core tends to retain an atmospheric mass of hydrogen com- 
parable to its own mass (MacLeod et al. 2012). Additionally, 
giant stars do not frequently get deposited into highly-bound 
orbits in which the core itself is likely to be disrupted, as 
the densities of their cores are ^10^ times larger than their 
envelopes, and their orbital migration into the loss-cone is 
largely dictated by diffusion (Wang & Merritt 2004; MacLeod 
etal. 2012). 

We can thus conclude that while helium-rich disruptions 
will occur occasionally, they will not be the dominant con- 
tributor to the rate, and as a result, it is highly unlikely by 
chance that these disruptions would be among the first to be 
observed. 

7.2. Inclusion of priors 

In this work, we do not make any assumptions about the 
distribution of any of our input parameters, with the excep- 
tion of our intrinsic variability parameter ay (see Section 5.1). 
For example, the distribution of stars around SMBHs is likely 
to possess a current mass function (CMF) that is strongly re- 
lated to the initial mass function (IMF), which would sug- 
gest that the most likely stars to be disrupted are those with 
M* ~ O.IM0 (Kroupa 2001). Additionally, we might expect 
that that grazing encounters (e.g. small P) should outnum- 
ber deep encounters (Frank & Rees 1976), and the black hole 
mass should follow established Mh-L relations (Graham & 
Scott 2013). 

This prior information could be used to further constrain the 
parameters of any given event. However, each of these priors 
has a great deal of uncertainty associated with it. In our own 
galactic center, it is not clear if the distribution of stars is sim- 
ilar to the general IMF observed in the field, especially given 
the prevalence of short-lived B -stars within several It-days of 
the black hole (Gillessen et al. 2009). These stars may have 
been deposited through binary disruption (Ginsburg & Loeb 
2006), or through a disk (Madigan et al. 201 1), both of which 
would lead to different distribution in /3 than what would 
be produced by a steady-state, spherically symmetric cluster 
around the black hole (a la Wang & Merritt 2004). Lastly, 
while a clear trend has been demonstrated between the lumi- 
nosity of the host galaxy and the mass of its black hole, there 
is significant scatter about this trend (Gultekin et al. 2009). In 
principle, once a significant number of TDEs have been iden- 
tified, these distributions can be determined from the collec- 
tion of fits to all disruptions, which could potentially improve 
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the accuracy of parameter estimations of future events. 

7.3. How BLRs can help us understand TDEs 

The (non-)existence of various Hnes in spectra acquired of 
TDEs can be used to great effect to constrain the parameter 
space of allowed encounters for any particular event. These 
features enable one to place a time-dependent size constraint 
on the size of the debris structure resulting from a tidal dis- 
ruption, which is directly related to the combination of three 
parameters: Mh, M*, and /3 (Equation 9). In this paper, we 
have focused specifically on Hell A4686 and Ha in regards 
to PSl-lOjh, but our technique could be used in general with 
other emission lines. PSl-lOjh appears to originate from a 
moderately-massive SMBH, but disruptions of stars by more 
or less massive black holes would respectively produce larger 
or smaller structures from which emission lines could be pro- 
duced. As an example, the disruption of the same star by a 
10^ Mq may show H/3 and Hel A5876 features early, with Ha 
appearing later, whereas a disruption by a lO^M© black hole 
may never show any helium or hydrogen emission lines. 

In the optical at z = 0, the number of strong emission lines 
is limited, but many more emission lines are available in the 
UV and X-ray where metals with larger ionization potentials 
begin to lose their electrons. These lines, which would be 
produced nearer to the SMBH, could potentially be used to 
constrain the size scale at early times, and would provide a 
spatial map of the accretion disk as it grows. It is critical that 
TDEs are identified early and followed up spectroscopically 
to obtain this valuable information. 

7.4. How TDEs can help us understand BLRs 

The BLR has long been used to measure the masses of black 
holes from the lag times observed in the response of various 
emission lines, which are thought to lie at various distances. 
However, there remains much uncertainty in these models, 
namely the form of the BLR itself. If the dissipation mech- 
anism within the debris disks resulting from tidal disruptions 
is similar to the dissipation mechanism that controls angular 
momentum transport in steadily- accreting AGN, it is reason- 
able to expect that the two structures should have many simi- 
larities in terms of their density and temperature profiles, ve- 
locity structures, and in the components of the structure that 
conspire to produce the emergent light. 

In this paper, we have made direct comparisons to BLRs 
in order to understand the emission features that are observed 
in a particular event. As we have shown, the dependence be- 
tween the distance at which a particular emission line is pro- 
duced and the flux originating from the central engine is even 
more exaggerated than in the case of steadily-accreting AGN, 
as some line-emitting regions do not exist at all due to the 
absence of mass beyond a certain distance. With a larger cat- 
alogue of TDEs, we can reverse the arguments presented here 
to learn more about the structure of the BLR present in TDE 
debris disks. 

7.5. Caveats and future directions 

At the time of this writing, PSl-lOjh is the only event that 
is claimed to be a TDE and also captures the rise, peak, and 
decay of the flare. By capturing all three phases, and with 
the addition of spectroscopic information, this event provides 
significantly more information on the underlying mechanisms 
than the small number of poorly sampled UV/optical TDEs 
that only capture the decay phase and may have no spectro- 
scopic data (Gezari et al. 2006, 2008; Cappelluti et al. 2009; 



van Velzen et al. 2011; Cenko et al. 2012). While the models 
presented here provide compelling evidence of the similarities 
between steadily-accreting AGN and luminous flares result- 
ing from the tidal disruptions of stars, there are many aspects 
that can be improved upon. Some uncertainties in the gen- 
eralized model presented here, such as details on the viscous 
processes that govern accretion and how matter light is repro- 
cessed, could potentially be resolved with a more-complete 
collection of well- sampled TDEs. 

It is clear from our purely hydrodynamical simulations that 
mere gas dynamics is incapable of generating the necessary 
dissipation for high mass-ratio encounters, as we described 
in Section 3.3. This suggests that magnetohydrodynamical 
simulations that focus on the nozzle region need to be per- 
formed to examine the growth of the MRI, which by our es- 
timate may be capable of providing the required dissipation. 
If this mechanism is incapable of operating, then it is possible 
that only deeply-penetrating encounters in which r^ ~ r^ will 
yield rapidly-rising light curves. 

A second critical uncertainty is our treatment of the repro- 
cessing layer, which is inextricably linked to the BLR of TDE 
debris disks. In this work, we have presumed that this repro- 
cessing layer is spherical, parameterized the amount of light 
absorbed by an average gray opacity, and have ignored po- 
tentially complex radiation transport and line-of- sight effects. 
It is also unlikely that the BLR relations we compare to here 
are identical for debris disks resulting from disruption, given 
their elliptical geometry and different radial mass distribu- 
tions. While the scaling relations determined for steadily- 
accreting AGN are likely to be similar to TDE scaling rela- 
tions, meaningful constraints on individual events can only be 
obtained by revising these relations to account for the differ- 
ences. 

Given a more-accurate prescription of how the viscous and 
reprocessing mechanisms operate, TDEFit can easily be im- 
proved to include these additional aspects of the problem, 
which can potentially yield very accurate estimates of the pa- 
rameters associated with individual disruption events. With 
a large library of TDEs, which will likely exist in the LSST 
era when potentially thousands of TDEs may be detected (van 
Velzen et al. 201 1), it should be possible to obtain detailed de- 
mographics of the stellar clusters that surround SMBHs. 

7.6. Lessons Learned 

For the readers convenience, we summarize the main find- 
ings of this paper below. 

1. The unbound material, while ejected at high veloc- 
ity from pericenter after a disruption, is gravitationally 
confined in the two directions transverse to its motion. 
This constricts the debris to a thin stream that presents a 
negligible surface area as compared to the emitting sur- 
face generated by the return of the stream to pericenter, 
and is unlikely to affect the flare's appearance. 

2. When material returns to pericenter, it is heated via hy- 
drodynamical shocks, but this dissipation is likely in- 
sufficient to explain the tight relationship between L 
and M for large-^ encounters. Additional dissipation 
via an MRI-like mechanism or through general rela- 
tivistic precession is probably required to explain this 
observed relationship. 

3. A disk that is truncated at 2rt fails drastically in explain- 
ing the observed flare, and cannot match the observed 
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shape of the Hght curves without extreme color evolu- 
tion. 

4. The light curve of PSl-lOjh is well modeled by a single 
blackbody whose temperature evolves weakly in time, 
and whose size is tens of times larger than r^. We spec- 
ulate that this distance is roughly co-spatial with the 
distance at which helium is doubly-ionized. 

5. The fact that Hell emission lines are observed, but Ha 
and H/3 are not, is consistent with the size constraint on 
the bound debris that is ejected from the nozzle point 
upon returning to pericenter. In general, the presence 
or absence of various emission lines can be used as a 
probe of the size of the elliptical debris disk. 

6. The parameters for PSl-lOjh of our highest- likelihood 
fits indicate that a O.56M0 main- sequence star was par- 
tially disrupted by a 1.1 x 10^ Mq black hole. As this 
involves the disruption of a common star by a common 
SMBH with an impact parameter near the expected av- 
erage, TDEs of the kind we associate with PSl-lOjh are 
likely to be among the most common sub-Eddington 
disruption events. However, given that we are analyz- 
ing a single event in this paper, we cannot eliminate the 



possibility that an event of this type was among the first 
observed due to observational bias. Once more well- 
sampled TDEs are available, a joint analysis of many 
events similar to what we perform here is required for 
a complete understanding of the demographics of tidal 
disruption. 
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